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1.  Introduction 


Over  the  past  60  years,  a  large  number  of  papers  have  been  written  about  the  properties  of 
graphene  (1 ),  due  not  only  to  its  defining  characteristic  as  the  ultimate  two-dimensional  (2-D) 
electronic  system  but  also  to  its  unique  band  structure  and  the  strongly  covalent  features  of  its 
atomic  bonding,  properties  that  it  inherits  from  the  semimetal  graphite.  Indeed,  the  earliest 
papers  on  graphene  (2)  treat  the  material  as  an  idealized  version  of  graphite,  i.e.,  a  perfect 
monolayer  of  carbon  with  a  hexagonal  graphite  structure.  Despite  the  apparent  simplicity  of  this 
“single-layer  graphite”  model,  the  picture  of  graphene  that  evolved  from  these  papers  is  one  of 
extreme  complexity.  This  is  because  the  relative  simplicity  of  its  constituents — atoms  of  a  light- 
atomic-number  element,  carbon,  in  a  simple  high-symmetry  lattice — favors  the  application  of  a 
high-powered  many-body  analysis  of  the  electronic  degrees  of  freedom,  which,  in  turn,  predicts 
ferromagnetism,  superconductivity,  charge-density  waves,  quantum  Hall  effects,  and  other  even 
more  novel  phenomena.  However,  in  view  of  the  extreme  difficulty  of  creating  such  monolayer 
structures,  many  of  the  theoretical  results  obtained  from  this  extensive  research  have  been 
considered  irrelevant  by  technologists.  As  has  often  happened  in  the  past,  experimental  advances 
have  now  overcome  the  fabrication  difficulties  to  the  point  that  the  reality  of  graphene  as  an 
object  of  laboratory  study  is  beyond  question,  and  so  these  theoretical  predictions  can  be  put  to 
the  test. 

On  the  other  hand,  the  technological  implications  of  graphene’s  physics — its  monatomic  nature, 
excellent  electron/hole  transport  properties,  and  mechanical  strength  among  others — have  driven 
an  effort  to  create  devices  that  (it  is  claimed)  could  rival  silicon-based  transistors,  and  possibly 
harness  the  various  exotic  properties  of  the  material  that  have  been  foreseen  by  theorists.  At  this 
time,  there  are  still  fonnidable  obstacles  to  making  graphene  devices  and  circuits  on  the  scale  of 
those  based  on  any  standard  semiconductor-based  material,  let  alone  highly  evolved  silicon 
devices.  However,  interest  in  graphene  devices  has  now  reached  a  point  where  models  of  such 
devices  are  needed  that  describe  the  terminal  characteristics  of  circuit  components  without  the 
need  to  track  the  interior  details  of  the  device  function.  The  availability  of  fast  computer 
simulations  of  ordinary  semiconductor  devices  has  motivated  efforts  to  incorporate  such  details 
in  the  device  computations,  but  the  problem  of  long  running  times  for  codes  based  on  these 
simulations  has  not  been  solved. 

In  this  report,  I  argue  that,  rather  than  imitating  those  models  of  standard  semiconductor-based 
devices  that  incorporate  massive  ab  initio  computer  calculations — metal-oxide  semiconductor 
field-effect  transistors  (MOSFETs),  high  electron  mobility  transistors  (HEMTs),  metal 
semiconductor  field-effect  transistors  (MESFETs),  and  their  ilk — the  development  of  graphene 
device  models  should  be  evolutionary,  i.e.,  the  models  should  first  mimic  the  simplest  standard 
models.  As  a  first  step  in  this  program,  I  discuss  the  simplest  semiconductor  device — the  field- 
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effect  transistor  (FET) — treated  at  the  lowest  level  of  device  modeling,  i.e.,  Shockley’s  gradual- 
channel  model.  My  choice  of  FETs  is  driven  by  ongoing  interest  in  radio  frequency  (RF)  devices 
here  at  the  U.S.  Army  Research  Laboratory  (ARL)  and  the  mature  understanding  of  the  high- 
frequency  behavior  of  FETs  in  general.  In  the  course  of  this  analysis,  I  highlight  places  where  the 
unique  properties  of  graphene  affect  the  predicted  terminal  characteristics  and  S-parameters  of  a 
possible  graphene  FET. 


2.  Description  of  Gradual-Channel  Modeling 


This  type  of  device  modeling,  which  was  developed  by  Shockley  (3)  in  the  pre-computer  era,  is 
based  on  highly  idealized  pictures  of  the  interior  fields  and  currents  of  the  device  under  study. 
For  a  FET,  the  standard  geometry  is  shown  in  figure  1.  The  device  consists  of  an  insulating 
substrate  for  electrical  isolation  and  mechanical  support;  a  conducting  “active”  layer  of  width  h 
whose  conductance  can  be  modulated  by  varying  the  2-D  charge  density  and  drift  velocity  of 
carriers  in  the  layer;  “source”  and  “drain”  electrodes  that  contact  the  active  layer;  a  dielectric  or 
depleted  non-conducting  layer  of  width  d  on  top  of  the  active  layer;  and  a  metallic  gate  electrode 
on  top  of  this  last  layer.  Like  any  transistor,  the  FET  has  three  terminals,  but  the  most  common 
configuration  for  its  use  is  as  a  circuit  two-port,  i.e.,  two  of  the  electrodes  (usually  the  source  and 
gate)  are  connected  together.  In  this  configuration,  the  FET  converts  an  input  current-voltage 
(I-V)  pair  (gate  voltage,  gate  current)  into  an  output  I-V  pair  (drain  voltage,  drain  current). 
Hence,  the  modeler’s  task  is  to  determine  two  equations  for  the  output  electrical  quantities, 
referred  to  as  “terminal  relations.” 


Figure  1.  Typical  FET  geometry  (gallium  arsenide  [GaAs]  MESFET). 
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When  voltages  are  applied  to  the  electrodes,  a  highly  non-unifonn  electric  field  exists  inside  the 
dielectric  layer  and  a  current  flows  within  the  active  layer.  The  device  operation  is  predicated  on 
the  degree  of  control  each  electrode  voltage  has  over  the  electric  field  and  current:  the  charge 
density  in  the  channel  is  enhanced  or  suppressed  by  that  “portion”  of  the  electric  field  controlled 
by  the  gate  electrode,  while  the  drain  electrode  controls  the  carrier  velocities  through  its 
“portion.”  While  the  concept  of  a  “portion”  is  necessarily  vague,  it  can  be  made  more  precise  by 
introducing  mutual  capacitances  (active  and  passive)  between  the  electrodes.  Because  the 
electrodes  are  asymmetrically  placed,  it  is  expected  that  the  current,  like  the  field,  is  non- 
uniform.  The  distribution  of  potential  within  the  structure  is  schematically  shown  in  figure  2. 


Figure  2.  Relation  between  internal  potential  and  terminal  voltages.  The  red  line  indicates  the  FET  channel,  the 
blue  line  the  gate. 

In  figure  2,  V (x)  denotes  the  voltage  at  point  x  along  the  graphene  layer  and  Vg,  V^,  and 

denote  the  voltages  on  the  source,  gate,  and  drain  electrodes,  respectively. 

The  guiding  physical  principle  behind  gradual-channel  FET  modeling  is  that  the  current  flowing 
from  source  to  drain  be  one-dimensional.  Charge  conservation  then  forces  this  current  to  be 
uniform  everywhere.  When  this  is  so,  the  constant  source-drain  current  I ^  Is  given  by 

ID=ZhJ  =  Zqhp(V)v(E)  (1) 

where  Z  is  the  width  of  the  FET  channel,  h  is  its  thickness,  and  the  (fluid-mechanical)  carrier 
current  density  J  =  qp(V)  v(E ) .  In  this  expression,  E  is  the  electric  field  along  the  channel,  q  is 

the  electron  charge,  0(C)  is  the  drift  velocity  of  the  carriers  in  the  field  E,  p  (V)  is  the  three- 

dimensional  (3-D)  number  density  (algebraic)  of  charge  carriers  in  the  channel,  and  V  is  the 
potential  “associated”  with  E.  For  extremely  thin  channels,  it  is  correct  to  introduce  a  2-D 
density  of  charge  carriers  in  the  channel  n  (F)  =  hp  (F) ,  so  that  the  channel  current  of  the  FET  is 

I D  =Zqv(E)n(V). 

It  is  clear  that  the  field  E  and  potential  V  are  defined  in  a  markedly  inconsistent  way.  For  the 
current  to  be  constant,  there  must  be  no  forces  perpendicular  to  channel,  implying  that  the  field  E 
must  point  only  from  source  to  drain  inside  the  active  layer.  Flowever,  if  this  E  is  the  gradient  of 


3 


some  potential  V,  the  latter  cannot  depend  on  the  voltage  applied  to  the  gate  electrode,  which 
contradicts  the  ability  of  the  gate  to  control  the  channel  charge.  This  difficulty  is  avoided  in  the 
gradual-channel  model  by  assuming  that  the  internal  electric  field  changes  direction  abruptly  at 
the  boundary  between  the  active  layer  and  the  dielectric/depleted  layer.  The  interior  field/current 
configuration  is  most  easily  represented  by  the  circuit  diagram  shown  in  figure  3,  which  is  that  of 
a  lossy  and  non-unifonn  transmission  line. 


Figure  3.  Lossy  transmission  line  model  of  FET  channel.  Resistors  are  functions  of  the  local  potential  difference 
between  channel  and  gate. 


Although  this  picture  of  the  internal  electric  field  and  carrier  transport  in  a  FET  is  fundamentally 
unrealistic,  it  has  proved  quite  successful  in  reproducing  the  FET  terminal  characteristics,  which 
points  to  considerable  insensitivity  of  these  characteristics  to  the  detailed  internal  device 
dynamics.  For  this  reason,  most  of  the  device  models  described  in  Sze’s  book  on  semiconductor 
devices  ( 4 )  start  with  this  model,  which  makes  it  an  excellent  template  for  comparing  graphene 
devices  with  their  better-known  counterparts. 


3.  Modeling  Methodology 


Once  the  general  concepts  of  the  gradual-channel  approach  are  accepted,  the  physics  of  the 
device  operation  are  needed,  that  is,  we  must  specify  the  functions  n(V)  and  u(is)  .  The  former 

is  derived  primarily  from  the  geometric  structure  of  the  device,  the  latter  from  microscopic 
features  of  the  carrier  transport. 

3.1  Channel  Carrier  Density 

Since  V  varies  along  the  channel,  modeling  n(V)  involves  the  dielectric  properties  of  the  device 

interior.  The  capacitor-transmission  line  picture  suggests  that  for  a  dielectric  layer  between  gate 
and  active  layer, 
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where  d  is  the  dielectric  layer  thickness  and  V  ^  is  a  threshold  voltage  below  which  there  are  no 

carriers  in  the  channel.  In  the  MOSFET  literature,  this  is  simply  known  as  the  “linear  model.”  Its 
use  implies  that  the  tenninal  characteristics  exhibit  no  nonlinear  behavior  with  drain  voltage,  and 
indeed  the  MOSFET  at  this  level  is  simply  a  gate-controlled  resistor.  In  an  effort  to  explain 
current  saturation  in  MOSFETs,  Hofstein  and  Heiman  (5)  added  a  variation  of  n  ( V )  with  the 

channel  voltage  V: 


Vr  -  V ,  -V 
d  \  G  th 


'(n-^vc-v^y 


(3) 


Within  the  gradual  channel  approximation  this  assumption  leads  to  the  I-V  curve 
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(4) 


Here  we  note  the  appearance  of  a  rather  embarrassing  feature  of  the  gradual-channel  model:  for 
values  of  VD  greater  than  V(  .  -  V ^ ,  the  current  I ^  decreases  and  eventually  changes  sign  at 

2{Vq-  Vth  j  -  Because  such  a  decrease  in  current  with  drain  voltage  under  DC  conditions  is  a  sign 

of  thermodynamic  instability,  its  appearance  here  must  be  unphysical.  It  is  tempting  to  assume 
that  it  arises  solely  from  the  poor  description  of  the  internal  field  in  the  gradual  channel  model, 
which  should  be  fully  2-D,  i.e.,  a  solution  of  the  2-D  Poisson  equation;  however,  at  least  one  of 
the  many  gradual  channel  models  for  MOSFETs — the  Pao-Sah  model  ( 6 ) — avoids  this  problem 
completely  without  introducing  2-D  effects.  In  this  model,  the  functional  form  of  n  (V )  is  quite 

complicated:  for  a p  —  MOSFET  with  semiconductor  3-D  bulk  charge  densities  p  ^  (majority) 
and  n ^  (minority), 


n(V)  =  —n  j  — - - 

V  2  P°  DV|/  F( pv|/,  V,A) 


(5) 


F{ Pv,  V,A)  = 
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U  P  1  A  P0  r 
where  P  =  -  — ,  A  = - ,  L  = 

kBT  PpO  ° 
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is  the  Debye  length  for  the  bulk  material, 


=  p  ^is  the  majority-carrier  bulk  charge  density,  and  i)/^  is  the  surface  potential,  which  is 
evaluated  from  the  gate  voltage  equation: 
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Since  its  introduction  in  1966,  this  model  has  regarded  by  designers  of  silicon-based  device  as 
too  awkward  and  time-consuming  to  use.  However,  from  the  standpoint  of  designing  graphene- 
based  devices,  its  formulation  includes  two  interesting  features:  in-channel  ambipolarity,  i.e., 
majority  and  minority  carriers  switch  places  within  a  pinched-off  channel;  and  diffusion,  which 
may  be  the  key  to  its  success  in  avoiding  the  need  for  fully  2-D  computations.  Unfortunately,  the 
treatment  of  carriers  in  the  channel  is  based  on  the  Boltzmann  equation  for  a  gas  of  classical 
charged  particles,  which  makes  it  unsuitable  for  the  analysis  of  a  truly  2-D  channel,  e.g.,  that  of  a 
HEMT. 


3.2  Drift  Velocity 

At  low  electric  fields,  the  drift  velocity  v(E)  of  current  carriers  in  a  semiconductor  can  easily  be 
derived  from  Ohm’s  law,  i.e.,  v(E)  =  p E  ,  where  the  mobility  p  is  taken  to  be  a  material  constant. 

However,  because  most  useful  devices,  including  FETs,  must  operate  under  large-signal 
conditions,  this  linear  result  must  be  generalized  to  include  large  fields.  One  way  to  do  this  is  to 
write  u  ( £ )  =  p  ( £ )  £ ,  where  p ( E )  is  derived  from  transport  theory  in  bulk  material.  Because  the 

derivation  of  p(E')  is  a  difficult  problem  even  in  this  case,  it  is  common  practice  to  specify  it 

empirically  (7).  From  the  standpoint  of  device  design,  the  most  important  requirement  for  such 
an  empirical  description  is  that  p(£)  incorporate  the  experimentally  observed  saturation  of  the 

carrier  drift  velocity  o (£)  in  most  semiconductor  materials.  For  graphene  FETs,  this  exercise  is 

highly  useful.  Several  types  of  drift  velocity  relations  are  considered  here,  among  them 


3.2.1  The  “Simulator”  Model 
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This  expression  (7)  is  commonly  used  for  silicon  (Si)  devices  and  was  in  fact  used  recently  by 
Meric  et  al.  (5)  to  model  graphene.  In  appendix  A,  it  is  shown  that  it  leads  to  the  I-V  relation 
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Using  the  MOSFET  linear  model  for  the  channel,  we  can  find  the  I-V  relation  explicitly. 
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Like  the  MOSFET  model  of  Hofstein  and  Heiman,  this  model  exhibits  the  problem  that  for 
values  of  VD  greater  than  y  the  current  I ^  decreases,  and  eventually  changes  sign  at  2V^ . 

This  shows  that  the  introduction  of  velocity  saturation  is  not  an  automatic  cure  for  the  problem  of 
thennodynamic  instability. 

The  expression  for  p(£)  given  previously  is  often  used  in  numerical  simulations,  but  its 

empirical  nature  suggests  that  is  has  very  little  physical  content.  Device  models  based  on  hot- 
electron  transport,  of  which  we  will  discuss  two,  are  much  more  interesting: 

3.2.2  Electron  Heating  Due  to  Phonons 

The  following  velocity-field  relation,  which  is  derived  from  inelastic  scattering  by  optical 
phonons  in  Conwell  (9),  is  occasionally  used  in  more  detailed  numerical  simulations: 
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A  derivation  of  the  I-V  curve  that  corresponds  to  this  expression  is  given  in  appendix  A  .  This 
expression  differs  markedly  from  the  “simulator”  result: 
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A  related  but  slightly  different  velocity-field  relation,  derived  from  inelastic  scattering  by 
acoustic  phonons: 
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Despite  the  similarity  of  equations  1 1  and  12,  the  latter  gives  a  rather  strange  result  for  the  I-V 
curve: 


These  two  expressions  will  be  discussed  in  a  separate  report. 

4.  Gradual-Channel  Models  for  Graphene 

With  regard  to  device  modeling  in  general,  there  are  four  peculiarities  of  graphene  that 
distinguish  it  from  bulk  materials: 

1 .  two-dimensionality 

2.  ambipolarity 

3.  zero  band  gap 

4.  a  conical  2-D  band  structure 

In  this  section  of  the  report,  I  examine  the  impact  of  these  features  on  three  kinds  of  FET  models, 
which  differ  in  the  way  velocity  saturation  is  modeled.  My  analysis  focuses  on  the  stability  of  the 
I-V  curves  (10),  which  is  an  important  criterion  for  the  usefulness  of  any  FET  model.  I  define  a 
DC  I-V  curve  as  stable  if  it  has  no  regions  where  drain  current  decreases  with  drain  voltage. 
Stated  another  way,  a  stable  DC  I-V  curve  should  exhibit  no  regions  with  negative  AC  output 
conductance.  Since  this  parameter,  along  with  the  transconductance,  are  needed  to  specify  the 
full  low-frequency  AC  response  of  the  FET  and  the  AC  equivalent  circuit,  stability  of  the  I-V 
curve  impacts  predictions  of  microwave  performance,  i.e.,  the  S-parameters  of  the  device. 
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4.1  Modeling  by  Meric  et  al. 


The  first  step  to  modeling  a  real  FET  structure  is  to  incorporate  a  varying  channel  voltage  into 
the  expression  for  the  channel  charge  density  n(V) .  In  their  paper  (5),  Meric  et  al.  provide  such 


an  expression  and  use  it  in  their  initial  effort  to  model  the  microwave  response  of  a  graphene 
FET.  Because  the  authors  wanted  to  study  the  role  of  ambipolarity  in  the  functioning  of  their 
device,  they  equipped  it  with  both  a  top  and  a  back  gate  (figure  4),  where  the  latter  is  used  to 
“dope”  the  FET  by  electrostatically  generating  either  n-  or  />typc  graphene  concentrations 
analogous  to  MOSFET  channels.  In  what  follows,  I  will  discuss  this  FET  model  using  data  from 
their  paper,  denoting  voltages  on  the  top  and  back  gates  by  Vgst  and  Vgsb,  respectively. 
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Figure  4.  Dual-gate  graphene  FET  structure  of  Meric  et  al. 


In  providing  an  expression  for  the  channel  carrier  density  n  (E) ,  Meric  et  al.  incorporate  an 

additional  peculiarity  exhibited  by  graphene  films:  the  existence  of  a  “minimum  conductivity”  in 
the  graphene  channel,  due  apparently  to  ambipolarity.  They  do  this  by  assuming  that  the 
minimum  conductivity  corresponds  to  a  minimum  carrier  density  n q  ,  and  then  using  the 

following  (empirical)  expression  for  the  function  «(E) ,  which  is  now  the  total  number  of  mobile 
carriers: 
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The  fact  that  this  expression  has  a  minimum  at  some  voltage  obscures  the  fact  that  neither  the 
number  of  electrons  nor  the  number  of  holes  has  such  a  minimum.  It  therefore  sidesteps  much  of 
the  discussion  of  ambipolar  transport,  but  this  is  legitimate  if  electrons  and  holes  have  identical 
mobilities,  which  may  be  the  case  in  graphene.  The  expression  contains  as  fitting  parameters  two 


threshold  gate  voltages  V  ® 


V. 
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gst 


at  which  the  minimum  carrier  density  n ^  is  reached  as  the 


back  and  top  gates  are  varied  separately.  The  design  numbers  for  this  FET  are 


Z  =  2.1  pm 
L  =  1.0  pm 


VF 

n0 


=  550 


cm 
V*s 
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=  1.x  10  cm/s 


=  5x10^  cm  2 


Cgsf  =552  nF/cm2 

VgtQ  =1.45  V 

C  ,  =  12  nF/cm2 
gsb 

V  ,°  =2.7  V 
gb 


(15) 


where  Z  is  the  top  gate  width,  L  its  length,  p()  the  channel  mobility,  u^,  the  Fenni  velocity  for 
an  infinite  graphene  sheet,  « ^  the  minimum  channel  density,  cgst  the  capacitance  between  the 
channel  and  top  gate,  and  c  /  the  capacitance  between  the  channel  and  back  gate.  Using  these 


numbers  and  a  back  gate  voltage  of  -40  V,  they  arrive  at  a  threshold  voltage  for  the  channel  of 
Vth  ~  2-38  V  ’ 


To  model  velocity  saturation,  they  use  the  simulator  expression,  which  leads  to  the  following 
formula  for  /  derived  in  appendix  A: 
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Although  they  have  written  several  papers  relating  to  this  modeling  effort,  Meric  et  al.  are 
unclear  as  to  the  value  of  the  saturation  velocity  i  >sa[ ,  although  they  mention  the  number 
7 

5.5  x  10  cm/s  (8).  They  justify  their  choice  of  the  simulator  expression  (equation  16)  implicitly 
by  referencing  to  reference  7,  but  presumably  simplicity  was  a  factor  as  well.  However,  because 
this  expression  was  developed  for  bulk-semiconductor  devices,  its  use  for  graphene  is  open  to 
criticism.  One  consequence  of  its  use  is  the  prediction  of  I-V  curves  that  have  a  point  of 
inflection,  or  “kink,”  when  the  biases  drive  the  Dirac  point  into  the  channel.  Their  experimental 
curves  do  indeed  show  such  a  kink;  however,  the  extent  to  which  the  latter  can  be  explained  by  a 
gradual-channel  model  needs  testing. 

Because  the  expression  for  the  I-V  curves  is  explicit,  it  is  easy  to  apply  a  stability  test  to  it  by 
computing  the  differential  output  conductance.  From  appendix  B,  we  find  that 
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where  I ^  is  given  by  equation  16.  In  appendix  B,  it  is  shown  that  the  condition  that  this  quantity 
always  be  positive  is  that 


Vds>0^Z‘>usal"(Vds)>/D 
Vds<0="  Zqvsatn(Vds)>  ~!D 


It  is  clear  that  if  this  inequality  is  violated,  the  violation  will  occur  within  a  certain  voltage 


interval  [F_  ,  V+  ]  bounded  by  points  where 


dl 
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=  0 ,  points  that  mark  a  peak  and 


r*  =  r± 


an  adjacent  trough  in  the  S-shaped  I-V  curve.  As  the  gate  voltage  is  varied,  this  interval  shrinks 
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to  zero,  at  which  point  V_=  V+.  Denoting  the  common  value  of  V_  and  V+  by  V ^  ,  we  see  that 
two  equations 
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are  satisfied  there  simultaneously,  so  that  the  curve  becomes  almost  flat.  The  second  equation  is 
the  condition  for  a  point  of  inflection  in  the  I-V  curve,  i.e.,  a  “kink.”  It  should  be  emphasized 
that  equation  19  and  equation  18  can  hold  at  the  same  time;  indeed,  these  conditions  are  required 
for  the  presence  of  a  physically  acceptable  kink.  If  F_  ^  V+  there  will  be  two  points  where 
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Let  us  demonstrate  this  for  the  expression  used  in  reference  8.  Since 
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we  find  explicitly  that 
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Figure  5  shows  plots  of  I  ^  versus  and  for  three  values  of  top-gate  voltage  for  the 

transistor  specified  by  the  numbers  given  in  reference  8  and  listed  above,  along  with  a  (rather 
7 

low)  value  of  2.5  x  10  cm/ s  for  u sa! .  For  this  transistor  the  threshold  top-gate  voltage  for 
negative  conductance  is  3.58  V.  The  blue  curves  are  below  the  threshold  voltage  for  negative 
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conductance,  and  show  a  kink  with  positive  slope;  the  red  curves  are  for  a  device  gated  at 
threshold,  where  the  kink  has  zero  slope;  and  the  yellow  curves  are  for  a  device  gated  above 
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are  shown  in  figure  6;  the  absolute  value  of  g  is  plotted  in 


gst 


accordance  with  reference  8. 


Figure  5.  (a)  ID  vs.  Vj,  for  the  simulator  model  and  (b)  /D  vs.  Vgst  for  the  simulator  model. 


For  the  real  devices  Meric  et  al.  investigated,  it  was  necessary  to  include  series  “access” 
resistances  associated  with  the  source  and  drain  contacts,  leading  to  the  following  modified  I-V 
expression  (11): 
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which  is  now  an  implicit  equation  for  the  current.  From  appendix  C,  we  find  that 
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Then  the  condition  that  this  quantity  always  be  positive  is  that 
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which  clearly  reduces  to  equation  18  for  no  access  resistances. 
Using  the  expression  introduced  in  reference  8: 
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we  again  find  explicitly  that 
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and 


V  =  si”h"1  M V  - + Vo  -  g* )  V = si'V  M  V  -  Vo  -  7* ) (29) 

Figure  7  again  shows  plots  of  1  versus  g,  and  g  for  the  same  three  values  of  top-gate 

D  ds  gst 

voltage,  for  the  same  transistor  and  numbers  given  above,  with  a  value  of  700  Q  for  the  access 
resistance  R^  .  Since  the  expression  for  l ^  involves  an  implicit  function  for  nonzero  R ^ , 

evaluation  of  1  requires  a  numerical  calculation;  once  it  is  computed,  however,  it  can  be 

d!  dl 

inserted  into  the  expressions  for — V_ .  The  output  conductance  s  ,  =  — —  and 

dV,  d  dV, 

ds 
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dl 

transconductance  g  = - are  shown  in  figure  8. 
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Figure  7.  (a)  ID  vs.  V*  for  the  simulator  model  with  access  resistances  and  (b)  ID  vs.  V*  for  the  simulator  model  with 
access  resistances. 
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Let  us  compare  the  results  of  this  calculation  with  the  results  for  no  access  resistance.  Note  that 
the  overall  current  output  is  suppressed,  which  is  unsurprising;  in  addition,  the  conductance  is 
much  larger,  a  more  obscure  consequence  of  the  additional  1400  Cl  of  access  resistance.  The 
threshold  top  gate  voltage  for  kink  formation  is  found  to  be  4.01  V,  an  increase  over  the  value  for 
no  access  resistance.  Figure  8(b)  shows  curves  similar  to  the  experimental  curves  given  in 
reference  12,  including  curvature  to  the  left  of  the  peak  that  is  absent  from  the  curves  for 
transistors  without  access  resistance. 


4.2  Corrected  Version  of  Simulator  Model 


The  authors  of  reference  12  mention  (and  then  seemingly  discard)  an  interesting  feature  of 
carrier  transport  in  graphene:  the  saturation  velocity  appears  to  be  a  function  of  carrier  density,  a 
phenomenon  that  has  never  been  reported  in  ordinary  FETs.  To  incorporate  this  effect,  they 
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parameters  (a=0.07).  The  use  of  this  function  in  simulations  is  problematic,  since  it  allows  o 


sat 


to  become  negative.  The  proper  form  of  the  I-V  curve  with  this  feature  included  is  derived  in 
appendix  B,  where  we  also  replace  the  semi-empirical  dependence  with  the  following  more 
reasonable  approximation:  since 


lim 
«— >0 


P 


r 


>£+!»  ^ 


we  can  replace  vsat  (ft)  =  »F  -j=-a  by  vsat  (ft)  = 

Vft 


P 


y 


1-^ 

v  P  j 

P°^7 

r( .  a 
\ln  l-l — \jn 


-  a , 


(30) 


with  little  error.  This 


V  P  J 

procedure  gives  us  a  well-behaved  dependence  on  n  everywhere  except  at  very  low  carrier 
densities. 
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The  derivation  given  in  appendix  A  is  correct  until  we  reach  equation  A-4.  In  order  to  include  the 
density  dependence  of  the  saturation  velocity,  we  need  the  expression 
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evaluated  in  the  previous  section: 
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where  0^v  =  sinh  1  >\vgst  ~  Vfh  )  ,  Qgd  =  sinh  1  k{vgst  ~  Vfh  -  ) .  The  second  integral  is 

more  challenging,  but  it  too  can  be  evaluated  in  closed  form: 
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is  the  inverse  function  of  the  Jacobi  elliptic  cosine  cn 


f  ,  o 

x  — 

V  2y 


For  the  stability  check,  we  take  the  same  derivative  as  before: 


dl D  Zqy.Q 


dV,  L 
as 


•(>*)■ 


D 

Zq$Up 


D 


n(Vds)  +  ^n(Vds)\  n(Vds)  Zqfiu 


1  +  - 


0 


0  a 


P*Jp£  [(3 


-h1+h2 


KY 


F 


l{Vds)+ln{Vds) 


0  a 


P vfl  [p 


-h1+h2 


®K) 


(36) 


If  we  calculate  the  I-V  curves  of  a  graphene  transistor  using  the  parameters  of  Meric  et  al.  and 
the  corrected  model,  we  get  the  curves  shown  in  figure  9. 


Figure  9.  (a)  ID  vs.  Vds  in  the  corrected  simulator  model  and  (b)  ID  vs.  in  the  corrected  simulator  model. 
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We  find  that  while  a  kink  is  still  present,  there  is  no  trace  of  negative  resistance  in  the  voltage 
intervals  shown.  This  feature  of  the  corrected  model  is  easier  to  see  from  a  detailed  comparison 
of  the  I-V  curves  shown  in  figures  10-16,  where  the  corrected  model  predicts  a  lower  current  at  a 
given  voltage  as  the  top-gate  bias  increases,  with  negligible  negative-slope  tendencies.  Further 
computations  indicate  that  the  onset  of  an  unphysical  region  in  the  I-V  curve  is  postponed  to  a 
higher  voltage  (5.21  V)  for  the  corrected  model,  but  the  actual  negative  slope  is  still  extremely 
small  for  voltages  as  high  as  7.58  V  (figure  8a).  This  indicates  that  the  corrected  version  is  a 
more  reliable  guide  to  the  transistor  behavior  over  the  range  of  drain  and  gate  voltages  of 
interest. 


Id,  mA 

o  < _ 

J  IT 

i  ^ 

l.U- 

n  5 

Vdc  \T 

0.0  0.5  1.0  1.5  2.0  2.5  3.0  ’ 

Figure  10.  IDv s.  Vds  for  V*  =  3.18  V,  blue  curve-uncorrected  model,  red 
curve-corrected  model. 


Figure  11.  Iq\ s.  Vds  for  V*  =  3.58  V,  blue  curve-uncorrected  model,  red 
curve-corrected  model. 
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Figure  12.  IDv s.  Vds  for  Vgsl  =  4.18  V,  blue  curve-uncorrected  model,  red 
curve-corrected  model 


Figure  13.  IDv s.  F*  for  V*  =  4.58  V,  blue  curve-uncorrected  model,  red 
curve-corrected  model. 


Figure  14.  ID  vs.  Vds  for  Vgst  =  5.18  V,  blue  curve-uncon'ected  model,  red 
curve-corrected  model. 
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Figure  15.  IDv s.  Vds  for  V*n  =  5.58  V,  blue  curve-uncorrected  model,  red 
curve-corrected  model. 


Figure  16.  IDv s.  Vds  for  Vgst=  7.58  V,  blue  curve-uncorrected  model,  red 
curve-corrected  model. 

Also  of  interest  are  the  output  conductance  and  transconductance  predicted  by  the  corrected 
model  compared  to  those  of  the  uncorrected  model,  which  are  shown  in  figures  17-22.  For  the 
transconductance  the  difference  between  the  two  models  is  not  large,  but  the  suppression  of 
negative  conductance  is  clearly  seen  in  the  output  conductance  for  the  corrected  model.  The 
more  rapid  decrease  in  the  conductance  away  from  the  origin  for  the  corrected  model  is  also 
apparent. 
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Figure  17.  g,„  vs.  Vds  for  Vgst=  3.18  V,  blue  curve-uncorrected  model,  red 
curve-corrected  model. 


Figure  18.  g,„  vs.  Vds  for  Vgs,  =  3.58  V,  blue  curve-uncorrected  model,  red 
curve-corrected  model. 


Figure  19.  gm  vs.  Vds  for  V*  =  4.18  V,  blue  curve-uncorrected  model,  red 
curve-corrected  model. 
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Figure  20.  sd  vs.  VD  for  Vgsl  =  3.18  V,  blue  curve-uncorrected  model,  red 
curve-corrected  model. 


Figure  21.  sD  vs.  Vds  for  Vs  3  .58  V,  blue  curve-uncorrected  model,  red 
curve-corrected  model. 
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Figure  22.  sD  vs.  Vds  for  V„v;  =  4. 18  V,  blue  curve-uncorrected  model,  red 
curve-corrected  model. 


5.  Conclusions 


It  is  interesting  that  the  simulator  model  used  by  Meric  et  al.  produces  a  negative  resistance  as 
soon  as  the  channel  current  exceeds  the  maximum  allowable  current  Zqvsatn(v ^  )  at  the  drain 

end  of  the  device.  In  standard  two-region  FET  models,  it  is  customary  when  this  current  is 
reached  in  part  of  the  channel  (the  “saturated  velocity  region”)  to  “pin”  the  variation  of  the 
channel  density  with  respect  to  beyond  this  point.  In  such  “two-region”  models  (13), 

increasing  then  increases  the  portion  of  the  channel  where  the  velocity  is  saturated,  which 

increases  the  current  slowly,  and  the  device  is  defined  to  be  “pinched  off’  once  the  entire 
channel  is  saturated  and  the  current  is  a  function  only  of  the  gate  voltage.  The  saturation  current 
that  flows  is  referred  to  as  1^  .  Two-region  models  do  not  exhibit  negative  resistance,  but  this  is 

because  of  the  ad-hoc  division  of  the  channel  into  segments. 

Thanks  to  the  papers  by  Meric  et  al.,  it  is  clear  that  simple  modeling  of  the  DC  and  microwave 
responses  of  graphene  FETs  is  feasible.  It  is,  however,  important  to  choose  a  model  that  includes 
as  many  of  the  peculiarities  of  graphene  as  possible,  in  order  to  avoid  a  misleading  picture  of 
device  function  leading  to  false  predictions  and  design  errors.  This  report  will  be  followed  by  a 
“physics-based”  modeling  effort,  in  which  the  model  described  here  will  be  “fleshed  out”  with 
more  realistic  features. 
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Appendix  A.  Derivation  of  I-V  Curves  From  Velocity-Field  Relations 


Three  types  of  drift  velocity  relations  are  considered  here 


A-l  The  “Simulator”  Form 


>(£)  = 


1  +  ^ 


(A-l) 


includes  an  absolute  value  function,  which  is  an  attempt  to  reflect  the  vector  character  of  the  drift 
velocity.  In  order  to  use  this  expression,  we  first  write  the  local  current  in  the  channel  at  point  x 


1 d  (x)  =  Zqv^E (x))n(V (x)) 


(A-2) 


Three  electrical  quantities  -  field,  potential,  and  current  -  are  connected  by  this  expression,  but  a 
second  relation  is  needed.  To  derive  it,  we  first  solve  equation  A-2  for  E.  The  absolute  value  in 
the  denominator  makes  this  slightly  tricky. 


,D=Zq — i+^ 

,,no  u 


Zq{i^E  n(V ) 


(A-3) 


Zq\x0n(V)\E\  =  !d 


since  the  denominator  in  equation  A-l  is  always  positive.  Then 
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Zq^n(y)\E\  =  \lD\  +  \lD 


ll0  |^| 


E  Z^(.Lq«(F)- 


^l0  lD 


=>  E  =■ 


z^0»(F)- 


^l0  ^Z) 


(A-4) 


From  this  we  get 


i+hs!£L,+ 

^  r/Y/ 


Z^h0»(F) 


^  7  tr/\  ^0  JD  7  (T,X  ^0  JD 

zW0n\V) - Zq\iQn{y) - 


J— =  1 _ M 

Pq  E  Zqn(V)\ 


(A- 5) 


We  now  integrate  equation  A-l  over  the  distance  x  from  source  to  drain,  whose  separation  we 
call  L: 


}  I ^dx  =  Zq\ 


j  ^  ^  j  ,  , 

0°  o,P'0l£l 


L  \JL  E  L 

1 — — TzrAv)dx  =  zd^\  1 

0i+V5  o[ 


Zqn(Vy 


i(V)Edx  (A-6) 


But  the  constancy  of  I ^  says  that  the  left  side  of  this  equation  is  just  I .  As  for  the  right  side, 


CIV 

note  that  E  = - =>  Edx  =  dV  ,  so  that 

dx 


no 

I  L  =  Zq\l  }  1- 


Zqn(V)  i 


i(V)dV 


where  is  the  source-drain  voltage. 


7  Vj  V, 

ds  ,  ,  Fn  “s  dV 

‘D~  f  n(V)dV  7 1 D  j  ZL 


L  '  u'  0  "so, 


v  v 

Zq\i„  ds  ds  dV 

Let  us  label  the  two  integrals  a  = -  j  n[V)dV  and  b  =  —  j  - 

L  0  L  0  vsat 

to  solve  is 


(A-7) 


so  that  the  equation 


ID=a~b,D 


(A-8) 
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The  integrals  a  and  b  have  positive  integrands;  hence,  they  have  the  same  sign  as  the  drain 
voltage  .  Figure  A-l  shows  plots  of  the  functions  f[ljj )  =  I ^  and  g[l ^  j  =  a  -b  I ^ 


Figure  A-l.  Plots  of  the  functions  /  (x)  =  x  and  g(x)  =  a -b\x\  for  (a)  a  =  3,  b  =  2  and  a  =  3,  b  =  0.4  and 
(b)  a  =  -3,  b=—  2  and  a  =  -3,  b  =  -0.4. 

For  Vj  >  0  (figure  A-l(a))  and  <  0  (figure  A-lb);  the  intersection  of  these  functions 

defines  the  solutions  to  equation  (A-3).  It  is  clear  from  figure  A- la  that  equation  A-8  always  has 
a  positive  solution,  and  in  addition  a  negative  solution  if  b  >  1;  likewise,  figure  A-lb  shows  that 
equation  A-8  always  has  a  negative  solution,  plus  a  positive  solution  if  b  >  1.  We  reject  the 
additional  solutions  as  unphysical,  in  which  case  we  easily  find  that 


'  D 


a 


1  +  b 

a 

T^b 


V,  >0 
as 


V,  <0 
as 


r  = 


d  i  +  \b\ 


and  so  the  I-V  curve  is 


Vds 

0  J  n(V)dV 


L 


!d  = 


0 


1  + 


v 

fi()  drS  dV 
L  0  vsat 


(A-9) 


(A- 10) 


In  a  bulk  semiconductor  the  saturation  velocity  is  a  material  quantity  and  therefore  constant  with 
voltage.  Then  it  comes  outside  the  integral,  and  so  we  have  simply 
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(A- 11) 


ID 


zq\- 10  vds 
—r -  I  n{V)dV 
L  0 


1  + 
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Zo 


sat 


V 


ds 


where  the  nonanalyticity  at  =  0  is  worth  noting. 

A-2  Hot-electron  Physics  of  the  Channel 

Hot-electron  physics  of  the  channel  leads  to  the  following  expression,  associated  with  scattering 
by  optical  phonons: 


o 


(£)=  0 


1  + 


f  M  £  ^ 

V 


(A- 12) 


This  resulting  I- V  curve  differs  markedly  from  the  form  derived  in  section  A- 1 .  Again  we  solve 
for  A: 
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and  as  before. 
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so  that 


r  = 
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Zqn(v)vSat 


n(V)dV 


(A- 15) 


A-3  Drift  Velocity  Associated  with  Acoustic  Phonon  Scattering 

For  the  drift  velocity  associated  with  acoustic  phonon  scattering: 
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(A- 16) 


u(£)  = 


m0e 


solving  for  E  gives 
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Since  Isat  >  I ^  ,  we  have 
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Appendix  B.  Stability  of  Meric’s  I-V  Curve  Without  Access  Resistance 


First,  we  calculate  the  output  conductance: 
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Then  for  F^  >  0  we  get  stability  when 
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while  for  V,  <  0  we  have 
ds 
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Appendix  C.  Small-Signal  Parameters  From  Meric’s  I-V  Curve  Without 
Access  Resistance 


C-l  Output  Conductance 
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C-2  Transconductance 

In  this  section,  we  make  the  dependence  of  n  on  both  gate  and  channel  voltages  explicit,  i.e.,  we 
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Appendix  D.  Stability  of  Meric’s  I-V  Curve  With  Access  Resistance 


Again,  we  calculate  the  output  conductance: 
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Appendix  E.  Small-Signal  Parameters  From  Meric’s  I-V  Curve  With  Access 
Resistance 
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Appendix  F.  Modeling  of  Graphene  FETS  Using  Meric  et  al.  With  Corrections 
(No  Access  Resistance) 


The  derivation  used  in  appendix  A  is  correct  until  we  reach  equation  A- 10.  For  tracking  the 
density  dependence  of  the  saturation  velocity,  we  need  the  following  expression: 
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Appendix  G.  Evaluation  of  Integrals  in  Electrical  Parameters  of  Graphene 
FETS  Using  The  Corrected  Model  of  Meric  et  al.  (No  Access  Resistance) 
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Since  n(V)  =  Jn  ^  +  (C^  /  qj  \V  -  Vq  j  ,  we  can  evaluate  the  integrals  that  appear  in  the 


FET  electrical  parameters: 
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Appendix  H.  Small  Signal  Parameters  for  Graphene  FETS  Using  the 
Corrected  Model  of  Meric  et  al.  (No  Access  Resistance) 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


2-D 

two-dimensional 

3-D 

three-dimensional 

AC 

alternating  current,  nonzero  frequency 

ARL 

U.S.  Army  Research  Laboratory 

DC 

direct  current,  zero-frequency 

FET 

field-effect  transistor 

HEMT 

high  electron  mobility  transistor 

I-V 

current  versus  voltage 

MESFET 

metal  semiconductor  field-effect  transistor 

MOSFET 

metal-oxide  semiconductor  field-effect  transistor 

RF 

radio  frequency 
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